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ABSTRACT 

Centaurus A, the closest active galaxy, has been detected from radio to 
high-energy gamma-rays. The synchrotron radiation by extremely high energy 
protons may be a suitable mechanism to explain the MeV to GeV emission 
detected by the instruments of the Compton Gamma-Ray Observatory, as 
coming from the inner jets. This scenario requires a relatively large magnetic 
field of about 10'' G that could be present only close to the central black 
hole. We investigate the spectral energy distribution (SED) resulting from 
a one-zone compact acceleration region, where both leptonic and hadronic 
relativistic populations arise. 

We present here results of such a model, where we have considered syn- 
chrotron radiation by primary electrons and protons, inverse Compton inter- 
actions, and gamma-ray emission originated by the inelastic hadronic interac- 
tions between relativistic protons and cold nuclei within the jets themselves. 
Photo- meson production by relativistic hadrons were also taken into account, 
as well as the effects of secondary particles injected by all interactions. The 
internal and external absorption of gamma rays is shown to be of great rel- 
evance to shape the observable SED, which was also recently constrained by 
the resuhs of HESS. 
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1 INTRODUCTION 

The elliptical galaxy NCG 5128 is the stellar body of the giant double radio source Centau- 
rus A (Cen A). The whole radio source extends ~ 10° on the sky, and it can be resolved 
down to sub-arcsecond scales in the inner radio structures which correspond to several pc 
because of its proximity (3.84 ± 0.35 Mpc, Rejkuba 2004). Cen A is one of the best exam- 
ples of a radio-loud AGN (a Fanar off- Riley Class I galaxy) viewed from the side (~ 70°) of 
the jet axis. The galaxy has an absorbing band of gas and dust projected across its stellar 
body. The nucleus ejects linear radio/X-ray jets, becoming sub-relativistic at a few parsec. 
At about 5 kpc from the core, the jets expand into plumes, and there are huge radio lobes 
that extend beyond the plumes out to 250 kpc. 

For further information, the reader is referred to the review by Israel (1998), which 

fl 

focuses on the observed properties of Cen A and their phenomenological interpretation □. 
Besides of the electromagnetic emission, Cen A has recently called the attention because of 
a striking clustering of ultra high-energy cosmic ray events observed by the Pierre Auger 
Observatory around its location (Abraham et al. 2007, 2008). This strongly suggests that 
protons can be efficiently accelerated up to very high energies in Cen A, as it has been 
suggested long time ago (Romero et al. 1996). 

Historically, Cen A has exhibited strong variability (more than one order of magnitude 
in flux) at X-rays (Bond et al. 1996). Noticeable variations are also present at radio wave- 
lengths, likely powered by accretion events. Very high resolution (VLBI) radio measurements 
are needed to separate the extremely compact nucleus from its surroundings. The difference 
between jet and counter-jet brightnesses can be explained by mildy relativistic Doppler 
beaming, which enhance the radiation of the approaching jet. From the observed brightness 
ratio, it appears that the northeastern jet is approaching and the southwestern jet receding 
at moderately relativistic speeds v ^ 0.45c (Jones et al. 1996; Tingay et al. 1998; see also 
Bao & Wiita 1997). 

Cen A has been observed repeatedly by all the instruments of the Compton Gamma-Ray 
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Observatory (CGRO). The spectrum determined by the Energetic Gamma- Ray Experiment 
Telescope (EGRET. Thompson et ah 1995) and the Compton Telescope (COMPTEL) can 
be fitted with a broken power law, with a spectral index steepening at high energies. The 
spectral energy distribution (SED) shows a broad peak around -Epeak ~ 0.1 MeV, with a lu- 
minosity of ~ 3 — 5 X 10^^ erg s~^ assuming a distance of 3.5 Mpc. At higher energies, only 
an integral flux upper limit, with a significance of ~ 0.4(T, was obtained by HESS through 
a short 4.2 hs live-time exposure (Aharonian et al. 2005). These observations, however, are 
not simultaneous, something that must be taken into account when pondering any model 
for the high-energy (HE) emission. 

Concerning the origin of this HE emission, leptonic models (e.g. an electron-positron 
beam as part of a two flow model) where the emission is dominated by Inverse Compton 
(IC), and synchrotron-self Compton mechanisms, have been suggested to explain the nuclear 
SED of Cen A. These models crucially require of low magnetic fields (< 10 G) to allow the 
electrons to reach high Lorentz factors (e.g. Marcowith et al. 1998, Ghisellini et al. 2005, 
Lenain et al. 2008). 

Several authors have proposed that injection of relativistic protons can take place in ac- 
tive galaxies (see Begelman, Rudak & Sikora 1990, and references therein). Models where the 
hadrons dominate the HE radiative outcome through photo-meson channels have been ap- 
plied to blazars (e.g. Mannhein et al. 1991, Mannheim 1993), and hybrid models for extended 
jet features in AGNs (Aharonian 2002, which includes the proton synchrotron component; 
see also Miicke, & Protheroe 2001, Reimer et al. 2004, and references therein). The radiating 
electrons can also be secondary particles produced in inelastic collisions by primary hadrons 
with ambient nuclei (e.g. Schuster et al. 2001). In that case neutrinos would be emitted 
along with the gamma rays. 

Here we have considered the emission of both hadronic and leptonic (primary and sec- 
ondary) particles in a proton dominated jet, similar to what has been applied recently to 
low-mass microquasars by Romero & Vila (2008). We aim to obtain information on the 
physical parameters related to the nuclear emission of Cen A. 
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2 OUTLINE OF THE SCENARIO 

We adopt an Eddington mass accretion rate onto the central supermassive black hole 
(M ~ 10*^ M©, Marconi et al. 2000). The accretion disc can be represented by the standard 
geometrically thin, optically thick, Shakura & Sunyaev model with a temperature profile 

Tdi.c(r) = 2.2 X lO^Ma/Vit^' K, (1) 

where the accretion rate is expressed in units of 10^^ g s~^, and the radius in 10^^ cm (Frank, 
King & Raine 2002). 



The inner edge of the disc is given by the last stable orbital radius (for a Schwarzschild 
black hole i?in = 6i?g, whereas for a Kerr black hole spinning in the same direction as the 
disc -Rin = Rg). The disc extends up to -Rout = 30/?^,. The gravitational radius of the black 
hole is Rg ~ 1.4 x 10^^ cm. The hot corona is modeled as a sphere with i?cor = emitting 
photons with spectral index a = 1.9, and with lower and higher cutoffs at 1 eV and 20 keV, 
respectively. Table [H presents the parameters of the model. Note that for the values assumed 
the Doppler factor is small 

^ 1. (2) 

The power redirected to the jets is a fraction of the accreted power, in accordance with 
the jet/disc symbiosis hypothesis: Lj ~ O.lMacc (e.g. Kording et al. 2006). About ~ 10% of 
the jet power goes to relativistic particles within a compact region close to the base of the 
jet, between zq = 50 Rg and Zf = 5zo (e.g. Bosch-Ramon et al. 2006, Romero & Vila 2008). 
The kinetic power of relativistic primary particles has two terms, corresponding to the both 
species originally injected (electrons and protons): U'^^ = U^^ + U^^. 

A proton dominated jet is quantified through the ratio a = U^^ / Uf^ ^ 100. The injection 
distribution follows 

Q,,p{E,z)<^E-h-^ (3) 

for the relativistic particles. Most of the jet content is in the form of a thermal plasma with 
a mildly relativistic bulk Lorentz factor F = 3. This plasma is roughly in equipartition with 
a tangled magnetic field Bq ~ 10^ G. The jet is assumed to expand in a conical way, thus 
the density of the cold material within it also decays as cx z~'^, and the energetic particles 
suffer the corresponding adiabatic energy losses. 
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Figure 1. Timescales for energy gain and losses of both protons (left) and electrons (right). The cooling rate by photo-meson 
production (p-f) includes only interactions with the proton synchrotron photons (the dominant field). 



We have obtained the steady state solution of the kinetic equations for the different 
species of particles taking into account adiabatic and radiative energy losses and the possi- 
ble escape of the particles from the acceleration region at a timescale tesc ~ (-2/ — -2o)/c- We 
neglected diffusion and convection effects in our treatment (see Khangulyan et al. 2007). 



Figure [2] illustrates the relevant timescales at the base of the jet. The synchrotron energy 
losses determine the cut-off in the energy distribution of the primary particles, E^^^. We 
obtain 

^max _ 2.2 X 10^^ eV, and = 4.7 x lO^eV. (4) 

At heights z > Zf the relativistic particles cool rapidly. 

In expression 1^, the injected relativistic particles have the canonic spectral index for 
standard first order Fermi diffusive acceleration. We have considered that such acceleration 
proceeds with an efficiency rj and on a timescale of tacc = E / [rjecBo) . If the steepening in the 
MeV-GeV spectrum of Cen A is considered as the cut-off signature of proton synchrotron 
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Table 1. Values for model parameters. 



Parameter: description [units] 


values 


• HIar'V }inlp masm 
ivi . Lfid^^iv injic uiOiSS livjy 1 






g 


-^out • outer accretion disc radius [Rg] 


30 


-^cor • radius of tlic corona [^?g] 


8 


zo'- jet initial point in the compact object RF [ijg] 


50 


z f : end of the acceleration region 


5zo 


B: viewing angle [°] 


70 


acor: corona photon index 


1.9 


X' jet semi-opening angle tangent 


0.1 


a: relativistic particles power-law index 


2 


V: macroscopic jet Lorentz factor 


3 


a: relativistic hadron-to-lepton ratio 


100 


L^"^: kinetic power in rel. particles [erg s~^] 


6 X 10*3 


7min: minimum Lorentz factor of rel. particles 


100 


r}: acceleration cfBciciicy 


10-4 


Bq: equipartition magnetic field [G] 


10-* 


no: number density of cold particles [cm~^] 


1.5 X 10^° 



Rg = GMjc^ 

radiation, since the magnetic field is fixed, we obtain a value r\ ~ 10~^ (see Aharonian 2000 
for details). The remaining free parameter, the minimum energy of the injected particles, is 
used to match the MeV luminosity inferred from the CGRO data. The normalization of the 
distribution is given by: 

gmax 

Lf = J d\ j dEEQ,{E,z), (5) 
Vol 

where the volume Vol is that of the acceleration region. 



3 THE HIGH-ENERGY EMISSION AND PROCESSING OF RADIATION 

The hadronic synchrotron radiation dominates the SED at MeV-GeV energies. The same 
primary relativistic protons produce pions through inelastic interactions with the cold mate- 
rial of the jet and with radiation fields. The target photons come from the accretion disc and 
the corona and are also produced in the jet itself mostly via the proton synchrotron radiation. 

Inelastic pp coUisions produce, through tt'^ decays, 7 rays in the 1 GeV- 10 TeV range. 
The computed luminosity is similar to the proton synchrotron one. The primary leptonic 
synchrotron and IC contributions result of minor relevance. 



The exact calculation of the interaction of the produced radiation fields with the relativis- 
tic particles and the radiation self-interaction is a complex problem. As emission is generated 
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in the same volume where it is strongly absorbed or reprocessed, a rigorous formalism should 
include the solution of a radiative transport equation. Concerning the high-energy photons, 
along the line of sight s, the specific intensity obeys 



where the photon emissivities qi come from the relevant radiative processes at energy E, 
that are usually not the same providing effective absorption, which occurs at energies of the 
order of em_in{E). At the same time, photon annihilation provides an injection term Qe that 
enters into the kinetic equation for the distribution of the leptons 



where dE/dt includes in-situ acceleration and energy losses (for the general case see, e.g., 
Ginzburg & Syrovatskii 1964). The cooling rates and emissivities of the IC and p7 pro- 
cesses are also determined by the properties of the radiation fields, therefore the differential 
equations ([6]) and ([7j) are coupled. Our current treatment is simplified by considering sub- 
populations that result from each term of the injection function, and the radiative emissiv- 
ities one-by-one. 

Whereas the decay of neutral pions leads to 7-ray emission, the decay of the charged 
pions originates a population of secondary electrons which extends to higher energies than 
the primary leptons (see Levinson 2006 and Orellana et al. 2007). The injected distribution 
function of these products can be computed using the expressions given by Kelner et al. 
(2006) and Kelner & Aharonian (2008). For the secondary electrons the kinetic equation to 
solve is slightly different from that of the primaries, resulting in a break of the distribution 
at E^^^ where the acceleration and cooling rates are balanced, and a high-energy cutoff that 
is a fraction of E^^^. 

In our model the proton synchrotron component extends down to optical wavelengths 
and provides a dense radiation field which suppress almost the whole pp contribution ex- 
cept for a tail of 7-ray photons with energies E > 10^^ eV. This internal energy-dependent 
absorption leads to the injection of relativistic secondary pairs with energies again higher 
than those of the primary leptons. Bottcher & Schlickeiser (1997) provide useful expressions 
to calculate the pair injection rates. Figure H] presents the obtained SED of the inner region 




(6) 




(7) 
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Figure 2. Computed optical depths, for internal absorption in the radiative field dominated by the p-synchrotron component, 
and for external fields taken as the sum of the terms due to the accretion disc and the corona. 



of Cen A. In the jets, the development of IC electromagnetic cascades is suppressed by the 
large magnetic field, because the secondaries are mainly cooled by synchrotron radiation 
(see, e.g., Aharonian, Khangulyan & Costamante, 2008). 

Outside the jet, in the vicinity of the compact 7-ray production region, the photon fields 
provided by the accretion disc and the corona absorb the emerging photons at energies 
greater than ~ 100 MeV. Figure [2] shows the dependence of the internal and external optical 
depths on the energy of the gamma rays. The large viewing angle of the jets (assumed to 
be perpendicular to the accretion disc) makes the path of the high-energy photons emerging 
from the jet to pass close to the disc. Hence geometrical considerations in the calculation 
of the optical depths can be relevant to shape the detectable spectrum. We have followed 
the treatment given by Becker & Kafatos (1995) to compute the opacity provided by the 
disc, and Dubus (2006) for the corona. Figure [3] shows the SED corrected by this external 
absorption plus the internal one. The data points corresponds just to the nuclear emission 
of the source (Lenain et a. 2008 and references therein). Observational constraints by HESS 
are also included but, as we mentioned, they were not simultaneously obtained. As Cen A 
is nearby source, the effect of absorption in the extragalactic background light can be safely 
neglected. In Figure H] we show the SED that results if only internal absorption is considered. 
Such a case might be relevant in case of advective inflows onto the black hole. 
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Figure 3. Computed spectral energy distribution for the inner region of Cen A. The grey band indicates the strong absorption 
by the dust lane, hiding the AGN core. The high-energy observational constraints are also shown. The data points corresponds 
to the core of the source. 



4 DISCUSSION 

The simplified model presented here is capable of reproducing the observed high-energy SED 
of Centaurus A as coming from a compact acceleration region at the base of its jets. The 
relativistic content is originally proton dominated but as secondary leptons are copiously 
produced and the model, as all so-called "hadronic" models, is actually lepto/hadronic. We 
have not considered here the effects of the magnetic field cooling of the decaying pions and 
muons during their short life-times, but the effects are not expected to be as dramatic as in 
the case of microquasars (see Reynoso & Romero 2008 for such a treatment). 



The low-energy component of the SED (not shown here) that is usually attributed to 
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Figure 4. SED emerging from the inner jets of Cen A, affected only by internal absorption. The contribution from different 
radiative processes considered are indicated. Sy stands for synchrotron radiation, whereas the subscript 1 or 2 indicates whether 
the particles are primaries or secondaries. 

the core of Cen A (Chiaberge et al. 2001) can result from the radiation of electrons injected 
in outer regions of the jets, where the magnetic field is lower. Such electrons plus "fresh" 
hadrons could be injected by the decay of neutrons which travel from the acceleration region 
without suffering synchrotron losses. Such neutrons are generated through the p + p — >■ 
p + n + T^'^ and p + 7 ^ n + vr"*" channels (e.g. Atoyan, 1992). The distance the neutrons 
travel before suffering /3-decay is given by 

dl3 = -fnCTn ^ (^J^^ PC (8) 

Thus the energy transport by neutral beams to the outer regions of the central source can 
lead to powerful contributions to the SED at low energies. The electrons injected from neu- 
tron decay cool through synchrotron radiation in regions where the field is now low. The 
electron population will be a power law that will mimic the neutron spectrum, which in 
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turn will reflect the original proton spectrum. The total power produce by these electrons 
from radio to optical will be ~ xq41-42 g-i^ accordance with the observations, for a 
magnetic field of < 1 G. Notice that the most energetic neutrons could reach distances of 
~ 100 pc. They will inject then energetic protons in the inner radio lobes, that could be 
accelerated then in successive steps up to very high energies (Romero et al. 1996). 

The physical parameters that were used for tunning our model to reproduce the CGRO 
observations are the acceleration efficiency rj and the minimum Lorentz factor of the rela- 
tivistic particles. The former resulted in a reasonable value of r/ ~ 10~^, compatible with 
acceleration in a parallel shock front (e.g. Protheroe 1998); and the latter is an order of 
magnitude lower than the one inferred in the SSC model of Chiaberge et al. (2001), namely 
7mm = 2 X 10^ 

Aharonian, Khangulyan & Costamante (2008) have pointed out that internal absorption 
could be a problem for the escape of high-energy 7-radiation from the production region. 
This absorption can lead to the formation of spectra with almost arbitrary slope in the 
TeV range. Our results follow that trend because the absorption model presented here has 
many free parameters, and in particular, those related to the absorbing external fields are 
essentially unknown since the core of Cen A is so well enclosed by the host galaxy. 

Following the literature we have considered a thin disc plus a hot corona to describe 
the emission of the accreting flow. But at the high accretion rate (Eddington) considered 
here, the actual case can be different, for example forming a thick disc as in Begelman & 
Meier (1982). Then advective solutions might exist for the inflow and the absorption could 
be much lower that what we estimated. 

Finally, if the strength of the magnetic field outside the jets is low enough, the devel- 
opment of IC electromagnetic cascades will be unavoidable given the great density of the 
radiation fields. Then the power contained in the photons with energy > 10^^ TeV that 
emerge from the jet can be reprocessed to lower energy photons, and the SED presented 
in Figure [3] should be taken as a lower estimate for energies i?^ > 1 GeV. Future, longer 
exposures with HESS or HESS II could detect the source, with a rather soft spectrum above 
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200 GeV. 
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